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We have calculated the atomic electric dipole moments (EDMs) d of ^He and ^^^Yb induced by 
their respective nuclear Schiff moments S. Our results are d('^He) = 8.3 x 10~^ and d('^^'^Yb) — —1.9 
in units 10~'^^S'/(e fm'') ecm. By considering the nuclear Schiff moments induced by the parity and 
time-reversal violating nucleon-nucleon interaction we find dC'^Yh) ~ O.GdC^^Rg). For =^He the 
nuclear EDM coupled with the hyperfine interaction gives a larger atomic EDM than the Schiff 
moment. The result for ^He is required for a neutron EDM experiment that is under development, 
where ''He is used as a comagnetometer. We find that the EDM for "'He is orders of magnitude 
smaller than the neutron EDM. The result for ^^^Yb is needed for the planning and interpretation 
of experiments that have been proposed to measure the EDM of this atom. 

PACS numbers: PACS: 32.80.Ys,31.15.Ar,21.10.Ky 



I. INTRODUCTION 

There are a number of experiments underway to 
measure the time-reversal violating electric dipole mo- 
ment (EDM) of the neutron and of various atoms and 
molecules. The measurement of a non-zero EDM would 
signal the presence of new sources of CP violation be- 
yond the standard model, standard model EDMs being 
undetectably small 

The best limit on the neutron EDM is \dn\ < 2.9 x 
10"^^ e cm (90% c.l.) [2] and that on an atomic EDM has 
been obtained for Hg, \dC^^llg)\ < 2.1 x lO-^^ecm (95% 
c.l.) 0]. The result for Hg constrains new physics scenar- 
ios largely in the nuclear sector, since Hg is a diamagnetic 
atom (total electronic angular momentum J = 0) and in 
lowest order the electric field couples to the nuclear spin. 
The Hg EDM is induced most efficiently from a nuclear 
Schiff moment; this is essentially a residual nuclear EDM 
(largely screened by atomic electrons) and is non-zero 
due to the finite size of the nucleus 4!\. 

In this work we calculate the atomic EDMs for ■^Hc 
and ^^^Yb induced by their respective nuclear Schiff mo- 
ments. The result for ^He is required for a neutron EDM 
experiment that is under development, where '^He is used 
as a comagnetometer 0] . There are proposals to mea- 
sure an atomic EDM in Yb 0, H, and preliminary 
work towards such a measurement is underway ^lOj . De- 
termination of the size of the EDM induced by the Schiff 
moment is needed to gauge the relative sensitivity of the 
experiment and will be required for interpretation of the 
measurement. We note that a calculation for the Yb 
EDM induced by the tensor parity and time-reversal vio- 
lating electron-nucleon interaction was performed in Ref. 
[m , however this is not a leading contribution in popular 
models of CP-violation 



II. METHOD OF CALCULATION 

The method we use in the current work is the same 
as one of the methods we used in our earlier work [131 



to calculate the EDMs of Hg, Xe, Rn, and Ra induced 
by nuclear Schiff moments. In the first method used in 
that work, which is the one adopted here, the atoms are 
considered as closed-shell atoms and the EDMs are calcu- 
lated using the relativistic Hartree-Fock (HF) and time- 
dependent Hartree-Fock (TDHF) methods. All atomic 
electrons contribute to the self-consistent potential. This 
approach is called the approximation, where is 
the total number of electrons. For the atoms Hg and Ra, 
another method was also used, where HF and TDHF cal- 
culations were performed for atomic cores consisting of 
N — 2 electrons, with the two valence electrons removed. 
The valence electrons were included at the next stage 
of the calculations where the configuration interaction, 
supplemented by many-body perturbation theory for va- 
lence electrons with the core, was implemented using the 
approach set out in [l3|. It was demonstrated in our 
work [l3] that the two methods give very close results 
for the respective EDMs of Hg and Ra. Therefore, in the 
present work we use only the simpler approximation. 
We briefly outline the method of calculation below. 

The atomic EDM induced in the many-body state N 
by the parity and time-reversal violating (P, T-odd) in- 
teraction HpT can be expressed as 
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where the sum M runs over a complete set of many-body 
states, En and Ej^i are atomic energies, and is the 
atomic electric dipole operator. The P, T-odd interaction 
Hamiltonian Hpx has the form 
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where ip is the electrostatic potential produced by the 
nuclear Schiff moment S which mixes states of opposite 
parity. The form for this potential that is suitable for 
relativistic atomic calculations is UM 
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where B = J p{R)R'^dR and p{R) is the nuclear density. 

In the V approximation, we can write the atomic 
EDM induced by the Schiff moment as 



d = 2^{SnpT\dz\n) , 



(4) 



where the sum runs over the relativistic Hartree-Fock 
core states \n), dz is the single-particle dipole operator, 
and \5npT) denotes the correction to the state \n) due to 
the P, T-odd Hamiltonian hpT- The correction \5npT) 
can be expressed as 



(5) 



where \a) corresponds to an excited state. It is found by 
solving the equation 



(/lo - tn)\5npT) = -hpT\n) 



(6) 



(Equivalently, one may calculate the correction to \n) 
from the electric dipole (El) field and take the matrix 
element of the weak Hamiltonian to obtain d.) 

In the approximation, polarization of the core due 
to the fields h pt and dz is accounted for by including the 
polarization due to one field using the TDHF method, 
e.g. by replacing hpx in (O by hpj, = hpx + SVpT, since 

J2n{^'^PT\dz\n) = J2ni^'>^PT\dz + SVd\n) . 

As a test of our wave functions we have performed 
calculations for the ionization potentials and the scalar 
polarizabilities a of the ground-state for each atom and 
compared them with the available experimental data. 
The polarizability has the same form as the EDM, EqlU 



-2E 



M 



\{N\Dz\M)\^' 
En — Em 



(7) 



with the operator HpT in ^ replaced by the dipole op- 
erator Dz- So the scalar polarizabilities are calculated 
simply by replacing the correction \ 6npT) due to the P, T- 
odd field by the correction \Snd) due the El field. 



III. RESULTS OF ATOMIC CALCULATIONS 

In Table U we present our results for the ionization po- 
tentials and scalar polarizabilities for He and Yb along- 
side the available experimental data. It is seen that the 
results for helium are in good agreement with experi- 
ment, while for Yb a deviation of about 15% is seen for 
the ionization potential. 

Results for the atomic EDMs for He and Yb induced 
by their respective nuclear Schiff moments are listed in 
Table im For comparison and easy reference we have also 
presented in the table the results from our previous work 
|12j . calculated in the approximation V'^ . The He EDM 
is very small, ^ 10~^^{S/{e fm"^)) e cm, and is stable 
with inclusion of core polarization, the value changing 
by ~ 10%. 



TABLE I: Ionization potentials (IP, in cm ^) 
larizabilities (q, in a%) for He and Yb 



and scalar po- 







IP 




a 






HF 


Exp. 


HF 


TDHF 


Exp. 


He 


201472 


198311" 


0.997 


1.32 


1.383" 


Yb 


43130 


50443'= 


124 


179 




"Ref. 












''Ref. 


m 










'^Ref. 


11 











As with Hg and Ra [12[, the largest contribution to 
the atomic EDM for Yb comes from the outer s elec- 
trons and this differs in sign to the overall contribution 
from the core. The final result for Yb is about 60% 
the size of the Hg EDM (in terms of their Schiff mo- 
ments). Our numerical result for Yb agrees with an es- 
timate obtained by scaling the Hg result, (i(^^^Yb) ~ 
d(i99Hg)(5Z2p)Yb/(5Z2i?)Hg ^ -1.65(Yb)/5(Hg), in 
the units of Table HIl Z is the nuclear charge, and R 
is a relativistic enhancement factor; see Refs. [1^ and 
[T^ [l9| , respectively, where the parametric dependence 
was found and more recently applied. 

One may be concerned about the huge corrections to 
Yb coming from inclusion of core polarization, the final 
value being over four times the HF value. We remind 
the reader that in our previous work [T^ we also saw 
such corrections (for Ra and to a lesser extent Hg, see 
Table HI]) . In Ref. [H] an entirely different method for 
the calculation of many-body corrections was also carried 
out and yielded results for the EDMs of Hg and Ra that 
differ from those obtained in the approximation by 
less than 10%. 

We expect that the TDHF value for the EDM of Yb is 
accurate to about 20-30%, the result for He more accu- 
rate. 



TABLE II: Atomic EDMs d induced by respective nuclear 
Schiff moments S in the HF and TDHF approximations, units 
are 10~^^(S'/(e fm^)) e cm. Results for He and Yb are from 
the current work; others are from our previous work Ref. [T^]. 



z 


Atom 


HF 


TDHF 


2 


He 


0.743 X 10"" 


0.826 X 10 


54 


Xe 


0.289 


0.378 


70 


Yb 


-0.416 


-1.91 


80 


Hg 


-1.19 


-2.97 


86 


Rn 


2.47 


3.33 


88 


Ra 


-1.85 


-8.23 



IV. ATOMIC EDMS IN TERMS OF THE 
P, T-ODD NUCLEON-NUCLEON INTERACTION 

We'd like to express the atomic EDMs for He and Yb 
in terms of a more fundamental parameter, in particu- 
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lar, the parameter specifying the strength of the P, T- 
violating nucleon-nucleon interaction rjj^j^i; we are inter- 
ested in this interaction because it leads to the largest 
Schiff moments. This will give a better idea of the rela- 
tive sensitivities of various atomic EDMs to fundamental 
physics. For instance, an order of magnitude enhance- 
ment of the nuclear Schiff moment may occur for de- 
formed nuclei with close levels of opposite parity [l8| 
(see also Refs. [1^ [2lj where nuclear EDM enhancement 
was considered); even more spectacular is the orders of 
magnitude enhancement that may arise in nuclei with 
octupole deformation [2^. 12^. 

The ground state of 7o^Yb has quantum numbers 
J'^ — 1/2^. While this nucleus is deformed, there are 
no opposite parity levels, 1/2+, close to the ground state 
and therefore there is no close- level enhancement. Nu- 
clei of ^^^Yb and ^^^Hg have the same quantum numbers 
and similar magnetic moments. In the shell model we 
consider the ground states to have an unpaired neutron 
in state p/i2- Therefore we may expect 

SC^^Yh) « SC^^Rg) = -1.4 X 10-^77„p e fm^ (8) 

so that the atomic EDM is 

d("iYb) « O.GdC^^Rg) « 3 X IQ-^'^Tjnp e cm . (9) 

The nuclear ground state for '^He is J"^ = 1/2+ and 
the magnetic moment is very close (~ 10%) to that 
of the neutron; we consider that in the ground state 
there is an unpaired neutron in the state si/2. Accord- 
ing to Ref. [l3|, the nuclear Schiff moment for '^He 
is S'(^He) ^ 0.1 X 10~®77„pecm. The Schiff moment 
scales with atomic mass A as A^^^. Scaling from ^^^Hg, 
Si^Ue) - S'(i99Hg) X (3/199)2/3 - 0.1 x lO-VpCcm, 
so the result of Ref. [l3| looks reasonable. The induced 
atomic EDM will then be d(^lle) « 1 x 10~30,y^p g qt[^ 

However, for very light atoms the finite size effect 
(Schiff moment) does not lead to the largest atomic 
EDMs. Another way to violate the screening of the nu- 
clear EDM is to take account of magnetic fields. Con- 
sidering the hyperfine interaction, an order of magnitude 
estimate for the atomic EDM arising from a nuclear EDM 
is H 

datom('He) - Za'^—dnncCHe) , (10) 



where Z is the nuclear charge, a is the fine structure 
constant, rrie and rup are electron and proton masses, 
and the nuclear EDM is denoted (inuc(^He). A more in- 
volved estimate was performed earlier by Schiff who ob- 
tained |d(3He)| = 1.5 X 10~^dnuc(^He) t^]. Using the 
value dnuc('^He) ^ 1 x lO^^^T^npCcm from Ref. [lM|, the 
size of the induced atomic EDM is then 

|d(^He)| - 1.5 X 10~2^r?„pecm, (11) 

larger than that induced by the nuclear Schiff moment. 

The size of the neutron EDM induced by the nucleon- 
nucleon interaction is dn = 0.5 x IO^^^t^^^ ecm [18|, and 
therefore it is seen that the atomic EDM of ^He is negli- 
gibly small compared to the neutron EDM d„, 

d{^Re) - 3 X IQ-^dri . (12) 



V. SUMMARY 

We have calculated the atomic EDMs for ^He and 
^^^Yb induced by their respective nuclear Schiff moments 
with the results (i(3He) = 8.3x 10"^ and d("iYb) = -1.9 
in units 10"^^(S'/(e fm'^)) e cm. The accuracy is about 
20-30% for Yb and is better for He. We also esti- 
mated the sizes of the nuclear Schiff moments induced by 
the P, T-violating nucleon-nucleon interaction. We find 
rf("iYb) - 0.6rf(i^^Hg). For ^He the nuclear EDM cou- 
pled with the hyperfine interaction gives a larger atomic 
EDM than the Schiff moment. Nevertheless, the helium 
EDM is orders of magnitude smaller than the neutron 
EDM and therefore may be neglected in the neutron ex- 
periment [1, Q . 
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